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Hamiltonian formulation of reduced Vlasov-Maxwell equations 
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The Hamiltonian formulation of the reduced Vlasov-Maxwell equations is expressed in terms of 
the macroscopic fields D and H. These macroscopic fields are themselves expressed in terms of the 
functional Lie-derivative £s generated by the functional <S, where £sJ~ = [S, J~] is expressed in 
terms of the Poisson bracket [ , ] for the exact Vlasov-Maxwell equations. Hence, the polarization 
vector P = (D — E)/47r and the magnetization vector M = (B H)/4-7r are defined in terms of the 
expressions 47rP = [S, E] + ■ ■ ■ and 47rM = — [5, B] + ■ ■ ■ , where lowest-order terms yield dipole 
contributions. 



I. INTRODUCTION 

Averaging dynamical equations over small spatial and 
temporal scales allows a simplified description and an 
easier numerical integration of these equations in order 
to identify the basic mechanisms at play. This averag- 
ing is traditionally performed by defining macroscopic 
quantities defined as integrals of the corresponding mi- 
croscopic quantities over small volumes and a small time 
interval. The macroscopic equations that are the dy- 
namical equations for the macroscopic quantities are also 
obtained by computing the integrals of the microscopic 
equations. Once obtained, these macroscopic equations 
are better suited to the investigation of the underlying 
processes at the space and time scales of interest. In gen- 
eral, this procedure suffers from three drawbacks: The 
first drawback is that the basic structure of the micro- 
scopic equations (e.g., Hamiltonian structure) might be 
lost in the integration process. The second drawback is 
that this procedure is associated with a loss of informa- 
tion, which means that once macroscopic quantities are 
computed (e.g., numerically), there is no way to go back 
to the microscopic ones. The third drawback involves the 
so-called closure problem: one may obtain a hierarchy of 
macroscopic equations that cannot be closed unless ad- 
ditional assumptions are included in the model. 

Canonical transformations allow one to perform such 
averages by offering a way to eliminate the dependency 
over fast variables. They have the significant advantage 
of not modifying the expression of the Poisson bracket. 
Another advantage is that the information on the mi- 
croscopic scales is encapsulated in the change of coordi- 
nates, and therefore can be recovered at any step in the 
reduction procedure, provided that the transformation is 
(locally) invertible. Such reduction procedures are par- 
ticularly useful in plasma physics. 

In magnetized plasma physics, where the high number 
of particles prevents a direct integration of the dynamics, 
the only hope is to rely on averaged equations. This is 
the purpose of guiding-center theory or gyrokinetic the- 
ory [ll 0] ■ As a consequence of the elimination of the fast 



time scales, associated to the rapid motion of particles 
around the magnetic field, in the Vlasov equation, polar- 
ization and magnetization effects arc introduced into the 
Maxwell equations. 

In this article, we propose a method for the direct 
Hamiltonian reduction of the Vlasov-Maxwell equations. 
The reduced equations are expressed in terms of the 
macroscopic fields 



D 
H 



E + 4?rP, 
B - 4ttM, 



(1) 

(2) 



which are obtained from the Vlasov-Maxwell Poisson 
bracket by the action of Lie transformations defined on 
a functional space. As an outcome of this procedure, we 
obtain expressions for the polarization P and magnetiza- 
tion M vectors. More precisely, we show how to define 
the fields D, H and F by a reduction process, so that 
Eqs. (|43| and (|45*|) are mapped into the reduced Maxwell 
equations 



<9D 



= cVxH - 47rJ 



dt 

VD = 4ttpr, 



R ■ 



(3) 
(4) 



where pr, and Jr are the reduced charge and current 
densities expressed in terms of moments of the reduced 
Vlasov distribution F. The remaining Maxwell equations 



dB 
~dt 
V-B 



cV X E, 







(5) 
(6) 



are left unchanged by the functional transformation. The 
main advantage of the Hamiltonian reduction is that 
we preserve the Hamiltonian structure of the Vlasov- 
Maxwell equations at each step of the process. 

In Sec. [H] we briefly review Lie transforms in the con- 
text of noncanonical and infinite dimensional Hamilto- 
nian systems. In Sec. IIIII we show that some of the re- 
sults obtained for particles are lifted to the Poisson alge- 
bra of functionals. In Sec. II VI we apply Lie transforms to 
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the Vlasov-Maxwell equations to derive the magnetiza- 
tion and polarization effects in magnetized plasmas. In 
particular we show that when the Lie transformation at 
the particle level does not depend on the electric field, 
the magnetic field is unchanged (no magnetization ef- 
fects) and the reduced distribution is identical to the one 
obtained in the case of external fields. 



II. FUNCTIONAL LIE TRANSFORMS 

The idea of dynamical reduction that we present here 
relies on the existence of a small parameter, e, and con- 
sists of identifying, for a given dynamical system, another 
system (which we refer to as reduced system) "close" to 
the original one, in the sense that the two systems become 
identical when e — 5-0. The reduced system possesses the 
desired properties, such as, for instance the existence of 
an invariant sub-space, so that the dynamics, up to some 
given order in e, is independent on a variable or on a 
parameter. 

Because our goal is to apply the method to Hamilto- 
nian systems, such as the Vlasov equation, a convenient 
mapping which permits to find the reduced dynamics, is 
the Lie transform. Indeed the latter possesses the use- 
ful properties of being invertible, close to the identity for 
e — > 0, and of being a Lie morphism p|. Because the 
systems of interest here are infinite-dimensional, the dy- 
namical variables being functions of coordinates z play- 
ing the role of labels, we make use of a Lie transform 
acting on Junctionals of the dynamical variables. Indeed, 
we consider a Poisson algebra consisting of functionals of 
some fields. We denote fa(z) the dynamical field vari- 
ables. The Poisson bracket is denoted [•,•], and the dy- 
namics of any functional J-[<fti(z), . . . , </>jv(z)] is given by 
the Hamiltonian H[fa(z), . . . , fa\(z)]: 



(7) 



The Poisson bracket satisfies the antisymmetry property 
[g, = - [J 7 , g], the Leibnitz rule [TQ, K] =T[G,K\ + 
[T, K] g, and the Jacobi identity 



t, [g, ic] + g, [k, t] + jc, [j?, g] = o, (8) 



where (J-,g,lC) are arbitrary functionals of the fields 

fa(z). 

We consider a given functional S of the fields fa, and 
consider the functional Lie transform generated by this 
functional acting on any other functional J- as 



~ e£s F = F-e[S, F] 



[S, [S, P[] 



The reduced functionals are defined by 



; J". 



(9) 



(10) 



Here there are two dynamics to be distinguished, even 
though they are generated using the same Poisson 



bracket [•,•]: the original dynamics generated by the 
Hamiltonian H and the reduced dynamics generated by 
the reduced Hamiltonian %. 

For the dynamical equations, we consider the reduced 
evolution operator defined by 



dt 



dt 



-e£* 



e e£s T, e e£s U 



From the property 

[e e£s T,e e£s g] = e e£s 



(11) 

(12) 



for any functionals T and g, we deduce that Eq. ([7]) 
becomes 



~~dt 



[J 7 , H] 



where the new Hamiltonian is defined by 



H 



dS 



s ds 



(13) 



+ .-., (14) 



with dS/dt viewed as a notation for [S,TL\. 

Since the present framework is that of noncanonical 
Hamiltonian systems, there exists in general a special 
class of Casimir invariants, which Poisson-commute with 
any observable J 7 . In other words, a Casimir invariant C 
satisfies [C, J 7 ] = for any functional T. These Casimir 
invariants are unchanged by functional Lie-transforms 
((TDD, i.e., C = e~ e£s C = C since [S,C] = 0. 

The reduced dynamics is linked to the dynamics of 
reduced field variables. First we define the reduced vari- 
ables : 



(15) 



Given the scalar invariance F{ipi) = J~(<fii), we recover 
Eq. (fTU)) when expressed in the original field variables fa . 
The dynamics of these reduced variables is the same as 
the reduced dynamics of the original variables, i.e., 



-e£s 



dfa _ defa 
dt ~ dt 



(16) 



where the left hand side of the equation contains the 
Lie transform in order to express the dynamics of the 
reduced field variables ipi into the old field variables fa. 
The above expression can be rewritten as 



-e£s 



[fa,H] e = [fa,H] 



(17) 



which is obvious according to the property (|12j) since 
[t/>i,it\ t — [e e£s fa , W] . Here the bracket [•, -] c denotes the 
Poisson bracket [•, •] expressed in the new field variables 
ifa- 
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III. CASE OF EXTERNAL 
ELECTROMAGNETIC FIELDS 

We consider the Vlasov equation in time-dependent 
external electromagnetic fields, for the particle density 
/(x, v,i) in phase space : 



| = _ v . V/ _i( E+ I xB )|. 
at m V c / ov 

The dynamics is generated by the Hamiltonian 



H(f,W,T) 



■fd e z 



w, 



(18) 



(19) 



where z = (x, v) and the extra variable w is due to the 
explicit time dependence of the fields (r is its canonically 
conjugate variable). The Poisson bracket between two 
functionals in the extended phase space (/; u;, r) is given 
by 



+ 



m j dv \ Sf dw dw Sf 
dT_dQ_ _ dT_dQ_ 
dw dr dr dw' 



(20) 



where the bracket {■, ■} is given by 



{/, 9} = 



1 



V/- 



B 



dg_ 
dv 
dj_ 

dv 



dv 
dg_ 

dv 



(21) 



that is to say that it is the noncanonical Poisson bracket 
for charged particles in a magnetic field, where the last 
term appears as a result of the use of the particle velocity 
v instead of the canonical momentum. 



A. Functional Lie-transform 

We perform a Lie transform generated by a functional 
S(f,w,T), so that the functionals J- are changed into 
J- = cxp(— e£$)F- The reduced variables (F,w,t) are 
transformed into 



F 



e e£s w, 



e e£s r. 



(22) 



We impose the natural condition that the time variable 
is unchanged (r = r), which implies that dS/dw = 0. 
The distribution function is mapped into 



F 



f 



f = f + e 
e [a, /} + ■ 



S, f 



(23) 



where the second expression shows the standard push- 
forward operation on the function /(z) generated by the 



function <j = SS/Sf and associated with the Poisson 
bracket {•,•}. 

For the reduction procedure, it is now more conve- 
nient to work at the level of the coordinates z of the 
particles. This is what is conventionally done in gyroki- 
netic theory |l|. More specifically, we present below an 
application of the proposed method to the reduction in 
the guiding center theory in the case of a static magnetic 
field and in the absence of electric field. In this example 
the reduction procedure consists in obtaining, under the 
assumption of high cyclotron frequency, a Vlasov system 
describing, up to a given order, the evolution of a guiding 
center distribution function, which represents the density 
in a phase space which is "smaller" than the initial one. 
More precisely, the reduction consists in eliminating the 
dependence on the particle gyration angle, which evolves 
on time scales much shorter than those of interest. The 
reduction procedure described in this example is based, 
to a great extent, on the analysis that was carried out 
by Littlejohn, at the level of particle dynamics in Refs. 

mm. 



B. Guiding-center Problem 

In a nutshell, Littlejohn [l5[ derived a sequence of 
(in general noncanonical) transformations that, starting 
from the canonical coordinates used to describe the mo- 
tion of a charged particle, led to a new set of coordinates, 
in terms of which, the corresponding dynamics is inde- 
pendent of the gyration angle, up to terms of the desired 
order in the small parameter e, which represents the ra- 
tio of the gyration period of the particle (with mass and 
charge set to unity) with respect to the characteristic 
time scale of the guiding center dynamics. For the sake 
of our example, we only consider the last step in the pro- 
cedure which performs the gyroangle average. We denote 
the new coordinates with (X, U, J, 9) and their relation 
with the original particle coordinates will be given below. 

The starting point for our example is the Hamiltonian 
functional (up to first order in e) 



H(F) = J d 6 ZF (i? + eHi 



(24) 



where H = \ U 2 + J B denotes the unperturbed Hamil- 
tonian, the first-order Hamiltonian is 

2 



(25) 



Hi = p ■ [U 1 k + - JVB 



1 



JU ( b • V X b — 2ai 



and the bracket (note that E = 0) 

[T,Q] = Jd 6 ZF |j>, g F } 



(26) 



with subscripts on functionals indicating functional 
derivatives. In Eq. (|2"5")) . p = (2J/B) 1//2 a denotes the 
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lowest-order gyroradius vector, k = b • Vb denotes the 
magnetic curvature, and a\ = —\ (c • Vb -a + a- Vb • c), 
with c defined such that a = b x c. The gyration angle 
of the particle motion is indicated with 9. An alternative 
orthonormal basis (1, 2, b) with a = 1 cos 9 — 2 sin 9, can 
be defined and through which the dependence of a and c 
on 9 can be made explicit. 

The distribution function F(X, U, J, 9) is the density 
in phase space identified by the coordinates (defined up 
to first order in e) 



X = x-ep 
U = u-e 



u k • Pq — Jq ( - b • V X b 



Oil 



J 



Jo 



1 



1 



D 



b-Vxb--p -VlnS 



(27) 
. (28) 

,(29) 
(30) 



where Jo = Vj_ /2B, x is the particle position, u = v • b 
and v± — |v x b| its parallel and perpendicular velocity, 
respectively, and ( is the physical (instantaneous) gyra- 
tion angle, and R = VI • 2 denotes the gyrogauge vector. 
The inner bracket in the bracket (l26l) is 



\Ji95gc — e \ 09 dJ at no 



d£dg 
dJ 89 



U 



B 



dX dU 



dg_dl 

ox ou 



b \dx ax 



(31) 



Because the magnetic field is external and no electric 
field is present, the observables of this Hamiltonian sys- 
tem consist of the functionals F(F), which depend on 
the distribution function F only. Notice that the Hamil- 
tonian (|24p is of the form where H\ depends explicitly on 
9. With the help of the functional Lie transform proce- 
dure described in Sec. lU we show how this dependence 
is removed at order 0(e) and pushed to terms of order 
0(e 2 ). 

According to Eq. (|10l) . the transformed Hamiltonian 
reads 

H(F) = e- e£s H(F) = H(F) - e[S(F), U{F)] + ■■■ 

= J d 6 Z(H + e(Hx - {S F , H } gc ) + ---)F 

(32) 



and comparing this last expression with (|32[) . one obtains 

H = H , (34) 
H 1 =H 1 -{S F ,H } gc . (35) 

Bearing in mind that we want Hi to be independent of 9, 
the first-order relation can be separated into ^-dependent 
and ^-independent parts: 



Hi = <#!>, 



= H 1 -{S F ,H o } i 



(36) 



(37) 



//, ./ T ( i b-V x b 



where (•) denotes ^-averaging and H = H — (H) denotes 
the ^-dependent part of H (it is assumed that (S F ) = 0). 
The condition (l36l) tells us that 



(38) 



On the other hand, the condition (1371) gives us an equa- 
tion to be solved with respect to S F - Its solution reads 

s ' = W-(T" + ! Jvh ' B ) + f °* (39) 

where ot\ = 80:2/ d9. Therefore, the functional that gen- 
erates the required Lie transform is 



S(F) = J d e ZFS F , 



with S F given by Eq. (|3U)) . The resulting Vlasov equation 
is then given by 



^ = [F,H(F)} = {H ,F} gc + e{Hi,F} 



(40) 



Because H and Hi are independent of 9, we identify 
an invariant sub-algebra, consisting of the functionals 
F((F)) of 9- independent distribution functions (F). 

Indeed, if the distribution function at t = is 9- 
independent, then, evolving according to the dynamics 
generated by [., H gc \, it will remain ^-independent at any 
time (since the Hamiltonian H. gc is also an element of 
this sub- algebra since it can be rewritten as a functional 
of (F)). This Lie sub- algebra realizes the dynamical re- 
duction. According to (fT6|) , the dynamics of the reduced 
distribution function F(X,U,J,9) = e e£s F(X 7 U 7 J,6), 
is given by 



9.1 



Of 



d e F 



(41) 



The generating function S is determined in such a way 
that it cancels all terms explicitly dependent on the gyro- 
angle 9 in Hi. Assuming that the transformed Hamilto- 
nian functional admits an expansion as 



H(F) = / d 6 Z(H Q + eHi)F + 0(e 2 ) 



(33) 



IV. REDUCED MAXWELL- VLASOV 
EQUATIONS 

The next example we consider is the case of a collision- 
less plasma whose dynamics is described by the Vlasov- 
Maxwell equations. These equations give the dynamics 
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of a distribution function /(x, v;t) of charged particles 
in phase space, together with an electric field E(x; t) and 
magnetic field B(x;t) : 

% = - W/^(e + ^xb)|, (42) 
at to V c / ov 



<9E 

at" 
<9B 

St" 



= cV X B - 4tt J, 

= - cV X E, 



(43) 
(44) 



where J = e J f\rd 3 v. The initial conditions of this dy- 
namical system are constrained by the two remaining 



Maxwell's equations 

V-E = 4ttp, 
V-B = 0, 

where p = e J fd 3 v, and a normalization condition 



d 6 zf = N. 



(45) 
(46) 



(47) 



The Vlasov-Maxwell equations (|421) -([44 j) are expressed 
in Hamiltoman form using the Vlasov-Maxwell Poisson 
bracket 



between two functionals T and Q of the Vlasov-Maxwell 
fields (/, E, B), and the Hamiltonian is 



H 



|v| 2 / d 6 z + 



fx 
"87 



(|Ef 



IBI 



(49) 



The first term in Eq. (|4"5)l involves the noncanonical Pois- 
son bracket (|2"Tj) . 

We perform a Lie transform generated by S(f, E, B) 
changing the functional T as T = exp( — e£s)J~- First 
we start by defining the reduced field variables D, H and 
F as 





(50) 



We notice the change of sign in the Lie transform which is 
a result of the scalar invariance J"(F, D, H) = F(f, E, B). 
Also we notice that D and H only depend on x and t since 
S has no explicit dependence on z = (x, v). 

We define the vectors P and M correspond to the po- 
larization and magnetization effects as 



(51) 





Equivalently, these vectors can be defined as 
P = 

M = — ( I 



i-(e^-l)E=i-[ < S,E] + 



1 



4tt 



)B 



1 



4tt 



We notice that V • M = since e e£s commutes with V. 
From p = V • E/ (4w), we have 



PR = V ■ 



E 

-in 



V- P 



-) 

4tt / 



P + V-P, 



which corresponds to the standard equation for the re- 
duced charge density (see Ref. @). The first-order ex- 
pressions for the polarization and magnetization are 



PlSi L [s , E ] = -cvxg 



Mi 



1 



[S,B] 



-cV x 



5S 



4tt L ' ' " " 8E 
where the displacement p\ is defined by 

i d rss 



J f Pl d 3 v, (52) 
(53) 



to dv \Sf 

which agrees with Ref. Q. The first-order expansion for 
F is given by 



Of J TO ov oE 



(54) 



Section Mil gives a natural explanation of the two terms 
in Eq. (|54|) : The first term of the right hand side of 
Eq. (|5"4"|) is associated with the action of the fields on the 
particles (like in the external field case), and the second 
one is intrinsic to the feedback from the particles to the 
fields. 

Next, we look at the reduced dynamics of the field 
variables (/, E, B) which is equivalent to looking at the 
dynamics of the reduced fields (F, D, H) as we have seen 
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in Sec. |TTJ First we consider the reduced dynamics of B 
and rewrite it as : 



~dt 



= [B,H] = c- e£s ([(e ££s - 1)B, H] + [B,H] 



In the first term of the rightmost hand side, we recognize 
the time derivative of M. By using Eqs. (|16[1 and ([51)1. 
we recover the dynamics of the reduced field H 

<9H „ „ A ( „ „ <9M\ 
— = -cV x D + 4?r cV x P - — — 
dt V dt J 

A similar derivation is done for the dynamics of the 
reduced field D by using the reduced dynamics for E 
rewritten as 



^ = [E, %] = e~ c£ * ([(e ejfi « - 1)E, U] + [E, H] 



Using Eqs. (fl6[) and ([oTp . the dynamics of D is given by 
<9D 



dt 



= cV x H - 47rJ ; 



where 



Jb = J - cV x M 



dP 

~dt~ 



In the case of external fields, a Lie transform in the 
functional space is associated with a phase space trans- 
formation at the particle level as it was seen in Sec. IIIII 



J 



Below we consider the case where the Lie transform is 
generated by S given by 



S(/,E,B) = j /(z)a(z;E(x),B(x))d 6 z, 



(55) 



such that SS/Sf = a. Furthermore we consider the case 
where the transformation done at the level of the par- 
ticles only depends on z and on B as it is done in the 
guiding-center theory. In this case we have <Sg = 0, and 
hence, from Eq. (|53p . the magnetization M is zero, which 
means that the magnetic field remains unchanged by the 
reduction procedure. In this particular case, we show 
that the distribution function is exactly the one obtained 
at the particle level by the Lie transform generated by a 
as in the case of external electromagnetic fields. First we 
notice that the right-most hand side of Eq. (j54j) vanished. 
Second, we notice that the Poisson bracket (l4"51) between 
two functionals T and Q of / and B gives 



Therefore we show by induction that 



F : 



7- 



where we recall that the Lie transformation is now gen- 
erated with the Poisson bracket {•,•} at the particle 
level. For instance the first order is already contained 
in Eq. ([54)) and the second-order term 



5, [SJ(z )]\ = y/(z)| ( 7(z;B(x)),^-y{/(z ), ( 7(zo;B(xo))}| ( i 6 



r 



Given that 
6 



{/(z ), <t(z ; B(x ))} = -{S b (z - z ), a(z; B(x))}, 



by linearity and since the function a does not depend 
on /. Using two integrations by part (noticing that 
J d 6 z{f,g} — for all function / and g of z), we ob- 
tain that 

S, [5,/(z )]] ={<7(zo;B( Xo )) ! {a(z ; B(x )), /(z )}} . 

When the Lie transformations at the particle level de- 
pend on the electric field, the passage from the particle 
description to the functional level, the generating func- 



tional is no longer as straightforward to construct from 
the generating function at the particle level. 
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